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Abstract 

We present the numbers of dimer- monomers on the Sierpinski gasket SGdin) at stage n with 
dimension d equal to two, three and four, and determine the asymptotic behaviors. The corre- 
sponding results on the generalized Sierpinski gasket SGd,b{n) with d = 2 and 6 = 3,4 are obtained. 

PACS numbers: 05.20.-y, 02.10.Ox 



* Electronic address: |scchang@mail.ncku.edu.tw| 
'''Electronic address: |lcchen@math.fju.edu.tw 



I. INTRODUCTION 



The enumeration of the number of dimer- monomers NohiiG) on a graph G is a classical 
model l|, y, 0]. In the model, each diatomic molecule is regarded as a dimer which occupies 
two adjacent sites of the graph. The sites that are not covered by any dimers are considered 
as occupied by monomers. Although the close-packed dimer problem on planar lattices has 
been expressed in closed- form almost half a century ago Q, [s], H] , the general dimer-monomer 
problem was shown to be computationally intractable [71]. Some recent studies on the enu- 
meration of close-packed dimer, single-monomer an d g eneral dimer-monomer problems on 
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It is of inter- 



regular lattices were carried out in Refs. j^, y, 1^, 
est to consider dimer-monomers on self-similar fractal lattices which have scaling invariance 
rather than translational invariance. Fractals are geometric structures of non-integer Haus- 
dorff d,„,e„sio„ real^ed bv repeated construction of an e.cnenta., shape on ptogtessive.y 



smaller length scales 19|, |20]. A well-known example of fractal is the Sierpins 



has been extensively studied in several contexts 
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We shall derive the recursion relations for the numbers of dimer-monomers on the Sier- 
pinski gasket with dimension equal to two, three and four, and determine the asymptotic 
growth constants. We shall also consider the number of dimer-monomers on the generalized 
Sierpinski gasket with dimension equal to two. 

II. PRELIMINARIES 



We first recall some relevant definitions in this section. A connected graph (without 



loops) G = (y,E) is defined by its vertex (site) and edge (bond) sets V and E (32|, |33 |. 
Let v{G) = \V\ be the number of vertices and e{G) = \E\ the number of edges in G. The 
degree or coordination number ki of a vertex Vi ^ V is the number of edges attached to it. 
A /c-regular graph is a graph with the property that each of its vertices has the same degree 
k. In general, one can associate monomer and dimer weights to each monomer and dimer 
connecting adjacent vertices (see, for example la]). For simplicity, all such weights are set 
to one throughout this paper. 

When the number of dimer-monomers N£,m{G) grows exponentially with v{G) as v{G) — *• 
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oo, there exists a constant zq describing this exponential growth: 

\tiNdm{G) 



zg 



hm 

t)(G)^oo 



v{G) 



(2.1) 



where G, when used as a subscript in this manner, imphcitly refers to the thermodynamic 
hmit. 

The construction of the two-dimensional Sierpinski gasket SG2{n) at stage n is shown 
in Fig. [H At stage n = 0, it is an equilateral triangle; while stage n + 1 is obtained by 
the juxtaposition of three ra-stage structures. In general, the Sierpinski gaskets SGd can be 



built in any Euclidean dimension d with fractal dimensionality D = ln{d+ l)/ln2 
the Sierpinski gasket SGd{n), the numbers of edges and vertices are given by 



22|. For 



eiSGdin)) 



(d + l) 



d 



{d+1 



\n+l 



v{SGd{n)) 



d+1 



[(d + l)" + l] 



(2.2) 



(2.3) 



Except the [d + 1) outmost vertices which have degree d, all other vertices of SGd{n) have 
degree 2d. In the large n limit, SGd is 2(i-regular. 
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SG2(0) 5G2(1) SG2{2) 




SG2{3) 



FIG. 1: The first four stages n = 0, 1, 2, 3 of the two-dimensional Sierpinski gasket SG2{n). 



The Sierpinski gasket can be generalized, denoted as SGd,b{n), by introducing the side 
length b which is an integer larger or equal to two [34] . The generalized Sierpinski gasket at 
stage n + 1 is constructed with b layers of stage n hypertetrahedrons. The two-dimensional 
SG2^b{n) with 6 = 3 at stage n = 1,2 and 6 = 4 at stage n = 1 are illustrated in Fig. [2l 
The ordinary Sierpinski gasket SGd{n) corresponds to the 6 = 2 case, where the index b is 
neglected for simplicity. The Hausdorff dimension for SGd,b is given by D = In / ln6 

34| . Notice that SGd,b is not fc-regular even in the thermodynamic limit. 
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FIG. 2: The generalized two-dimensional Sierpinski gasket SG2,b{n) with 6 = 3 at stage n — 1,2 and b — 4 
at stage n = 1. 

III. THE NUMBER OF DIMER-MONOMERS ON 5'G2(n) 

In this section we derive the asymptotic growth constant for the number of dimer- 
monomers on the two-dimensional Sierpinski gasket SG2{n) in detail. Let us start with 
the definitions of the quantities to be used. 

Definition III.l Consider the generalized two-dimensional Sierpinski gasket SG2,bin) at 
stage n. (a) Define M2fi{n) = N£)M{SG2,b{n)) as the number of dimer-monomers. (b) Define 
f2,b{n) as the number of dimer-monomers such that the three outmost vertices are occupied 
by monomers, (c) Define gr2fi{n), gl2,b{i^), gt2,b{n) as the numbers of dimer-monomers such 
that either rightmost, leftmost or topmost vertex, respectively, is occupied by a dimer and the 
other two outmost vertices are occupied by monomers, (d) Define hr2^b{n), hl2^b{n), ht2,b{n) 
as the numbers of dimer-monomers such that either rightmost, leftmost or topmost vertex, 
respectively, is occupied by a monomer and the other two outmost vertices are occupied by 
dimers. (e) Define t2fi{n) as the number of dimer-monomers such that all three outmost 
vertices are occupied by dimers. 

It is clear that the values gr2^b{n), gl2,b{i^), gt2,b{n) are the same because of rotation 
symmetry, and we define g2,hij') = gi^2,b{n) = ghfiiji) = gt2fi{n). Similarly, we define 
h2,b{n) = hr2fi{n) = hl2fi{n) = ht2,b{n). Since we only consider ordinary Sierpinski gasket in 
this section, we use the notations M2(n), f2{n), (72(^)5 h2{n), and t2{n) for simplicity. They 
are illustrated in Fig. [3], where only the outmost vertices are shown. It follows that 

M2(n) = /2(n) + 3g2{n) + 3/i2(n) + t2{n) (3.1) 
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for non-negative integer n. The initial values at stage zero are /2(0) = 1, ^'2(0) = 0, 
h2{0) = 1, ^2(0) = and M2(0) = 4. The values at stage one are /2(1) = 4, ^2(1) = 4, 
/i2(l) = 3, t2(l) = 2 and M2(l) = 27. The purpose of this section is to obtain the asymptotic 
behavior of M2(n) as follows. The five quantities M2(n), f2{'n), g2{'n), h2{n) and t2{n) satisfy 
recursion relations. 




M2{n) /2(n) gr2{n) ghin) 5*2 (n) 




hr2{n) hl2{n) ht2{n) t2{n) 



FIG. 3: Illustration for the configurations M2(n), f2{n), gr2{n), ghin), gt2{n), hr2{n), hl2{n), ht2{n) and 
t2{n). Only the three outmost vertices are shown explicitly for f2{n), g2{n), h2{n) and t2{n), where each 
open circle is occupied by a monomer and each solid circle is occupied by a dimer. 

Lemma III.l For any non-negative integer n, 

M2(n+1) = Ml{n)-3M2{n)[g2{n) + 2h2{n) + t2{n)]^ 

+3[h2{n) + t2{n)][g2{n) + 2h2{n) + ^2(^1)]^ - [h2{n) + t2{n)f , (3.2) 
f2{n + 1) = /|(n) + 6f^{n)g2{n) + 3f^{n)h2{n) + 9/2^52'^ + 252H 

+6f2{n)g2{n)h2{n) , (3.3) 

g2{n + 1) = f^{n)g2{n) + 2/|(n)/i2(n) + 4/2(n)^2'(^) + f^nMn) + 8/2(n)^2(n)/i2(n) 

+3glin) + 2f2{n)g2{n)t2{n) + 2f2{n)hl{n) + Agl{n)h2{n) , (3.4) 

/i2(n + 1) = f2{n)gl{n) + 4/2(71)^2 (n)/i2(n) + 2gl{n) + 2f2{n)g2{n)t2{n) + 7^^(n)/i2(n) 
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+Sf2in)hl{n) + 2f2{n)h2{n)t2{n) + 2gl{n)t2{n) + Ag2{n)hl{n) 



(3.5) 



t2(n + l) = c/^(n) + 6^^(n)/i2(n) + 3c/^(n)t2(^) + 9^2(ri)/i2('i) + 2/^2 




+Qg2{n)h2{n)t2{n) . 



(3.6) 



Proof The Sierpinski gaskets SG2{n + 1) is composed of three SG2{n) with three pairs 
of vertices identified. For the number M2{n + 1), the unallowable configurations are those 
with at least a pair of identified vertices originally occupied by dimers. Therefore, the three 
configuration with a pair of identified vertices occupied by dimers should be subtracted from 
all possible configurations M|(n). However, this procedure subtracts out also configurations 
with two pairs of identified vertices occupied by dimers that should be added back as illus- 
trated in Fig. m Finally, the configuration with three pairs of identified vertices occupied 
by dimers should be subtracted, and Eq. fl3.2p is verified. 



FIG. 4: Illustration for the expressfon of M2{n + 1). Certain outmost vertices of SG2{n) are not shown 
means they can be occupied by either dimers or monomers. 

As illustrated in Fig. [5], the number /2(r;, + 1) consists of (i) one configuration where all 
three of the SG2{n) are in the f2{n) status, (ii) six configurations where two of the SG2{n) 
are in the f2{n) status and the other one is in the g2{n) status, (iii) three configurations 
where two of the SG2{n) are in the f2{n) status and the other one is in the h2{n) status, 
(iv) nine configurations where one of the SG2{n) is in the f2{n) status and the other two 
are in the g2{n) status, (v) two configuration where all three of the SG2{n) are in the g2{'n) 
status, (vi) six configurations where one of the SG2{n) is in the f2{n) status, another one is 
in the g2{n) status and the other one is in the h2{n) status. Therefore, we have 



/2(n + 1) = /|(n) + 2fi{n)[gr2{n) + ghin) + gt2{n)] + fl{n)[hr2{n) + hhin) + ht2{n)] 
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+2f2{n)[gt2in)gr2{n) + gt2{n)gl2{n) + gr2{n)gl2{n)] 
+f2in)[grl{n) + gllin) + gtl{n)] + 2gr2{n)gl2{n)gt2{n) 
+f2{n)[gt2{n)hr2{n) + gl2{n)ht2{n) + gr2{n)hl2{n)] 

+f2{n)[ht2{n)gr2{n) + hl2{n)gt2{n) + hr2{n)gl2{n)] . (3.7) 

With the identity gr2{n) = ghiji) = gt2{n) = g2{n) and hr2{n) = hl2{n) = ht2{n) = h2{n), 
Eq. fl3.3p is verified. 




FIG. 5: Illustration for the expression of f2{n + 1). The multiplication of three on the right-hand-side 
corresponds to the three possible orientations of 56*2 (n +1). 

Similarly, gt2{n -|- 1) = g2{n + 1), ht2{n + 1) = h2{n + 1) and t2{n + 1) for SG2{n + 1) can 
be obtained with appropriate configurations of its three constituting SG2{n) as illustrated 
in Figs. O [7] and [8] to verify Eqs. (13.41) . (13. 5p and (13. 6p . respectively. 




FIG. 6: Illustration for the expression of gt2{n +1). The multiplication of two on the right-hand-side 
corresponds to the reflection symmetry with respect to the central vertical axis. 

Eq. (13. 2p can also be obtained by substituting Eqs. (13. 30 - 03. 6p into Eq. (13. ip . □ 
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FIG. 7: Illustration for the expression of ht2{n + 1). The multiplication of two on the right-hand-side 
corresponds to the reflection symmetry with respect to the central vertical axis. 



+ ^x3+ i-Ax3+ i-Ax3+ AAx3+ AAx3 










FIG. 8: Illustration for the expression of t2{n + 1). The multiplication of three on the right-hand-side 
corresponds to the three possible orientations of 56*2 + 1). 

There are always 27 = 3^ terms in Eqs. (13.31) . (13.41) . (13. 5p . (13.61) because there are three 
possible choices for each of the three pairs of identified vertices: both of them are originally 
occupied by monomers, or either one of them is originally occupied by a monomer while 
the other one by a dimer. The values of M2(n), f2{n), g2{n), h2{n), t2{n) for small n 
can be evaluated recursively by Eqs. (I3.2l) - (l3.6p as listed in Table [H These numbers grow 
exponentially, and do not have simple integer factorizations. To estimate the value of the 
asymptotic growth constant defined in Eq. (12.11) . we need the following lemmas. For the 
generalized two-dimensional Sierpinski gasket SG2,b{n), define the ratios 

J2A"^) 92,b{^) h2,b{n) 

For the ordinary Sierpinski gasket in this section ,they are simplified to be a2{n), (32{n) and 
72 H- 



Lemma III. 2 For any positive integer n, the magnitudes of f2{n), g2{n), h2{n), t2{n) are 
ordered as 

t2{n) < h^in) < g2{n) < h{n) , (3.9) 
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TABLE I: The first few values of M2{n), f2{n), g2{n), h2{n), t2{n). 



n 





1 


2 


3 


4 


M2{n) 


4 


27 


10,054 


499,058,851,840 


60,978,122,299,433,248,924,629,725,740,007,424 


f2{n) 


1 


4 


1,584 


78,721,368,064 


9,618,673,427,679,675,357,952,788,786,053,120 


92{n) 





4 


1,352 


66,974,056,448 


8,183,299,472,241,085,511,976,093,040,508,928 


/i2(n) 


1 


3 


1,148 


56,979,607,552 


6,962,123,286,110,084,944,276,569,997,705,216 


t2{n) 





2 


970 


48,476,491,776 


5,923,180,596,700,062,197,918,947,839,311,872 



then 



< i2{n) < I32{n) < a2{n) < 1 



(3.10) 



Proof Eq. fl3.9p is valid for the first few positive integer n by the numbers given in 
Tabled! By Eqs. (E3D-(ISSD, we have 



[a2{.ri) - I32{n)][2 + ha2{n) + 2a2(n)/32(n)] 



+ [1 + 2a2{n)][al{n) - I32{n)^2{n)\ 



(3.11) 



^2(^ + 1) h2{n + l) 



/|(n)c/2(n) /2(n)c/i(n) 



[a2{n) - I32{n)][2 + 3a2(n) + 3/52(ri) + 2a2{n)p2{n)] 



+I32{n)[a2{n) - 72(71) + 2[al{n) - /?2(n) 72(^1)]] , (3.12) 



h2{n + l) t2{n + l) 



/2(n)c/i(n) ^i(n) 



[«2(ri) - Hn)][2 + 6/32 (7i) + 2(3l{n) + 2/32(^)72(71)] 



+/32(n)[a2(n) - I2{n)][l + 2/^2 (n)] . 



(3.13) 



Eq. (13. 9p is proved by mathematical induction if Eqs. fl3.1ip -f l3TT3|) are larger or equal to 
zero. Equivalently, we need 72(77) < j32{n) < a2{n)^ which can be proved also by induction 
since the following two expressions are larger or equal to zero for any positive integer 77. 



92{n + 1) 
f^{n)g2{n) 



/2 ( 77 + 1) /72 ( 77 + 1) 
/|(77) /2 (77)^1(77) 



9 



= [a2{n) - p2{n)]^[l + Sa2{n) + ^^{n) + (n) [/^a (n) + 72(n)] + 4:al{n)pi{n)] 



+ [al{n) - Pi{n)][[a2{n) - 72(n)][l + 2a2(n)/32H] + 2al{n)[Pi{n) - jl{n)]] 



+ [a2{n) - P2{n)][4:a2{n)[al{n) - P2{n)j2{n)] + a2{n)P2{n)[a2{n) - j2{n)] 



+a2{n)p2{n)[P2{n) - 72 H]] + 4a2(^)[a2(^) - Pl{nW2in) - 72(71)] 



+ [alin) - /32(n)72(n)]2 + 2a2H[/?2H - 72(n)]2[2/?2(n) + aUn) + a2in)/3iin)] , 



(3.14) 



h2in + l) 
J2in)gl{n) 



92{n + l) t2{n + l) 
fi{n)g2{n) g'l{n) 



[Mn) - (32{n)Y[l + 4[32{n) + 272(71) + 9f3i{n) + Q[3i{n) + 6[3i{n)^2{n) + 2(3i{n)^^{n)] 



+ [a2{n) - P2{n)][2a2{n)[p2{n) - 72(72)] + 4pi{n)[a2{n) + p2{n) - 272(72)] 



+2a2in)Pi{n)[Piin) - ^Hn)] + 2(3l{n)[a2in)(32{n) - ^^.(n)]] 



+i32in)[alin) - Plin)][l32in) - 72(22)] [1 + 8/^2(22)] 



+P^in)[P2{n) - 72(n)]'[l + 2a2{n) + 2/^2(22) + Aal{n)] . 



(3.15) 



It is clear that < 72 (^) since all the quantities f2{n), g2{n), h2{n), t2{n) are positive. 
Finally, 02(22) < 1 once Eq. (13.91) is established. □ 



The values of 02(22), (32{n), l2{n) for small 22 are listed in Table HTl 



Lemma III. 3 Sequence {a2{n)}^^i decreases monotonically, while sequences {72 (t^)}^! 
and {t2(^)//2('^)}^i increase monotonically. The limits 02 = lini„_>oo ct2 ('^) ; P2 = 
liiRn^oo P2{n) , 72 = lim„_oo72('^) exist. 
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TABLE II: The first few values of 02 (n), /?2("')i 72("')- The last digits given are rounded off. 



n 


1 


2 


3 


4 


02 (n) 


1 


0.853535353535354 


0.850773533223540 


0.850772150002722 


Mn) 


0.75 


0.849112426035503 


0.850771337051088 


0.850772150002159 


72 (n) 


0.666666666666667 


0.844947735191638 


0.850769141078411 


0.850772150001597 



Proof By Eqs. (13. 3p and (13.41) . we have 

/2(n + l)g2in) - g2in + l)/2(n) 



/|(n)(72(n) [2/2(n) + 5^72^ + 2/^2 (n)] [asH - P2{n)] 



-[f2{n) + 2g2in)][ 



/IH g2in) 



> . 



(3.16) 



By induction, a2{n) decreases as positive n increases using the results of Lemma [III.2I By 
Eqs. (13. 5p and (13.61) . we have 

hin + l)/i2H - h2in + l)t2{n) 
= f2{n)h2in)[gj{n) + 4:g2{n)h2{n) + 3hl{n)][a2in) -72(n)] 
+2gl{n)h2{n)[g2{n) + 2h2{n) + t2{n)][(32{n) - 72(72)] 



+2f2{n)hl{n)[g2{n) + ^2H][||^ - |||] > , 
such that 72(^2) increases as positive n increases. Finally, we have 
t2in + l)f2{n)- f2{n + l)t2{n) 

= ^'^""^^W) ~ W)^ + Qf2{n)g2{n)h2{n)[a2{n) - 72H] 

+3fl{n)g2{n)t2{n)[a2{n) - P2{n)] + 9f2{n)gi{n)h2{n)[P2{n) - 72(n)] 



(3.17) 
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-2h{n)gl{n)Un)[ 



g2{n) hl{n) 



> 



(3.18) 



fiin) gl{n)t2{n) 

where the last inequahty holds because of the combination of the third and the last terms: 

hl{n) 



2f2{n)g2{n)t2{n) 



> 



(3.19) 



^f2{n) g2{n)' /2(n) gl{n)t2{n) 

Because the sequence a2{n) decreases monotonically and bounded below, the limit a2 
exists. Similarly, sequence 72 (^^) and t2{n)/f2{n) increases monotonically and bounded 
above so that the limits 72 and lim^^oo ^2(^)//2('^) exist. It follows that the limit 
lim^^oo ^2 in) I f2 (n) exists since 



lim 



t2(n) _ h2in) ^.^ t2{n) 



*oo f2{n) f2{n) n~*co h2{n) 

such that (3 = lim„_,oo ^2(^)/fl'2('^) exists because 



h2{n) 



hm 

n^oo f2[n) 



hm hm 



f2{n) 



92{n) 



□ 



(3.20) 



(3.21) 



With the existence of the limits 0^2, P2, 12, and 72 < /52 < 02, we have 

n^oo /2(n) 5(2(^+1) 

_ (1 + 3^2)^ + 2al + 3a2p2il + 2^2) 

l + 2p2 + 4^2 + P2I2 + 8a2p2 + 3ai + 2a2/9272 + 2a2/3| + 4:alp2 

by Eqs. fl3.3p and (13.41) . which leads to the following result. 



(3.22) 



Corollary III.l The three limits 72, (^2 and 02 are equal to each other. 



Lemma III. 4 The asymptotic growth constant for the number of dimer-monomers on 
SG2{n) is bounded: 

4_,„/,(„) + !W1±M!!!)1 < ,,,, < 4-,„/,(™) + Ml±^^ , (3.23) 

where m is a positive integer. 
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Proof Let us define A2(^) = f2{n + l)//|(n). By Eq. (13.30 . we have 

A2(n) = [1 + 3a2{n)]^ + 2al{n) + 3a2(n)/?2(n)[l + 2a2{n)] . (3.24) 
It is clear that 1 < A„, < 27, and 

[l + 272(m)]3 < [l + 272(n)]3 < [l + 2(32{n)f < X2{n) < [l + 2a2H]^ < [l + 2a2{m)f (3.25) 

for n > m. By Eqs. (12.31) and (13. ip . we have 

lnM2(n) ^ 2 ln[l + 3Q2(n) + 3a2{n)f32{n) + a2(T2)/32(n)72(T2)] 21n/2(n) 
^;(^G'2(n)) 3(3" + 1) 3(3" + 1)' ^ ' 

where 

ln/2(n) = lnA2(r2-l) + 31n/2(n-l) 

= lnA2(n- 1) + 31nA2(ri-2) + 3^1n/2(n-2) 



n-l 

= ^3"-i-^lnA2(j) + 3"-"^ln/2(m) (3.27) 

for any m < n. By the definition of the asymptotic growth constant in Eq. (12.11) . 

lnM2(n) 

ZSG2 = lim , , 

n^QO v[bG2{n)) 

— hm ^ ^"^^^ 3a2(^) + 3a;2(Ti)/92(ra) + a2{n) /32{n)'y2{n)] 

~ 3(3"+ 1) 

2 3"-i--'' In 

+ hm 



2 E;=™ 3"-i-^' In A2(j) + 2 [3"-- In Mm)] 



3(3" + 1) 

The proof is completed using the inequality (I3.25p . □ 

The difference between the upper and lower bounds for zsg2 quickly converges to zero as 
m increases, and we have the following proposition. 
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Proposition III.l The asymptotic growth constant for the number of dimer-monomers 
on the two-dimensional Sierpinski gasket SG2{n) in the large n limit is zsg2 = 
0.656294236916.... 

The numerical value of zsg2 can be calculated with more than a hundred significant 
figures accurate when m in Eq. (13.231) is equal to seven. 



IV. THE NUMBER OF DIMER-MONOMERS ON SG2fi{n) WITH & = 3,4 

The method given in the previous section can be applied to the number of dimer- 
monomers on SGd,b{f^) with larger values of d and h. The number of configurations to 
be considered increases as d and h increase, and the recursion relations must be derived 
individually for each d and h. In this section, we consider the generalized two-dimensional 
Sierpinski gasket SG2fi{n) with the number of layers b equal to three and four. For SG2,3{n), 
the numbers of edges and vertices are given by 

e(^G2,3(n)) = 3x6", (4.1) 
7 X 6" -I- R 

v{SG2An)) = ^^r^ , (4-2) 

5 

where the three outmost vertices have degree two. There are (6" — l)/5 vertices of SG2,3{n) 
with degree six and 6(6" — l)/5 vertices with degree four. By Definition IIII.ll the number 
of dimer-monomers is M2,3(n) = f2,3{n) +35'2,3('^) + 3/i2,3(^) +^2,3(^)- The initial values are 
the same as for SG2: /2,3(0) = 1, (72,3(0) = 0, /i2,3(0) = 1 and ^2,3(0) = 0. 

The recursion relations are lengthy and given in the appendix. Some values of M2_3(n), 
f2,3{'n), g2,3{n), h2^3{n), t2^3{n) are listed in Table UTTl These numbers grow exponentially, 
and do not have simple integer factorizations. 

The sequence of the ratio defined in Eq. (13.81) {«2,3('^)}$^i increases monotonically and 
{72,3('^)}J^Li decreases monotonically with < «2,3(?t.) < 72,3(?t-) < 1, in contrast to the 
results for SG2{n). The values of a2,3{n), P2,3{n), 72,3 (^) for small n are listed in Table ITVl 

By a similar argument as Lemma IIII.41 the asymptotic growth constant for the number 
of dimer-monomers on S'G'2,3(^) is bounded: 

51n/2,3M+6 ln[l+2a2,3M]+ln[l-+3a2,3M] ^ ^ ^ 5 In /2,3(m)+6 ln[l+272,3 (m)]+ln[l+372,3 (m)] / . o\ 



7x6" 
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TABLE III: The first few values of M2,3(n), f2,3in), g2,3{'n), /i2,3("')i t2,3in-)- 



n 





1 


2 


M2,3(n) 


4 


425 


755,290,432,490,932 


/2,3N 


1 


66 


116,464,644,336,176 


a2,3{n) 





56 


100,722,462,529,064 


h2,3{n) 


1 


49 


87,108,127,443,640 







44 


75,334,018,236,644 



TABLE IV: The first few values of a2,3("')5 /?2,3("')i 72,3 ('i-)- The last digits given are rounded off. 



n 


1 


2 


3 


Oi2,3{n) 


0.848484848484848 


0.864832955126948 


0.864833096846111 


/32,3(n) 


0.875 


0.864833178780796 


0.864833096846111 


72,3 (?^) 


0.897959183673469 


0.864833402432925 


0.864833096846111 



with m a positive integer. We have the following proposition. 

Proposition IV. 1 The asymptotic growth constant for the number of dimer-monomers on 
the generalized two-dimensional Sierpinski gasket SG2^3{n) in the large n limit is zsG2-i ~ 
0.671617161058.... 

The convergence of the upper and lower bounds remains quick. More than a hundred 
significant figures for zsg2,3 '^^^ ^e obtained when m in Eq. (14.31) is equal to five. 
For SG2,4:{n), the numbers of edges and vertices are given by 

e{SG2,4{n)) = 3 X 10" , (4.4) 
4 X 1 0" -I- 5 

v{SG2An)) = , (4.5) 

where again the three outmost vertices have degree two. There are (10" — l)/3 vertices of 
SG2,i{n) with degree six, and (10" — 1) vertices with degree four. By Definition IIILll the 
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number of dimer- monomers is M2,4(n) = f2,i{n) + 3(^2,4 (^) + 3/i2,4(^) + t2,4(^)- The initial 
values are the same as for SG2- /2,4(0) = 1, (72,4(0) = 0, /i2,4(0) = 1 and ^2,4(0) = 0. We write 
a computer program to obtain the recursion relations for SG2A{n). They are too lengthy 
to be included here and are available from the authors on request. Some values of M2,4(n), 
/2,4(^), g2,'i{n), /i2,4(n), t2,4('^) are listed in Table El These numbers grow exponentially, and 
do not have simple integer factorizations. 



TABLE V: The first few values of M2,4(n), /2,4(n), g2,4:{n), h2^i{n), t2^i{n). 



n 





1 


2 


M2,4(n) 


4 


14,278 


7,033,761,314,434,948,243,456,944,474,554,222,281,728 


/2,4W 


1 


2,220 


1,095,249,688,634,151,454,219,516,689,432,826,798,080 


52,4 (ra) 





1,914 


940,563,707,718,765,231,855,988,194,853,818,067,968 


h2,i{n) 


1 


1,640 


807,724,574,091,886,425,362,687,789,454,449,995,776 


t2,i{n) 





1,396 


693,646,780,368,841,817,581,399,832,196,591,292,416 



The sequence of the ratio defined in Eq. (13. 8p {02,4 (jt.)}^^^ decreases monotonically and 
{72,4(^)}J?Li increases monotonically with < 72,4('^) < a2,4(^) < 1, the same as the results 
for SG2{n). 

The values of a2,4(^), /32,4(^), 72,4('^) for small n are listed in Table IVTl 



TABLE VI: The first few values of 02,4(^)5 f^2,A{n)^ 72,4 ('t-)- The last digits given are rounded off. 



n 


1 


2 


3 


a2,4(ra) 


0.862162162162162 


0.858766468942539 


0.858766468941692 




0.856844305120167 


0.858766468941199 


0.858766468941692 




0.851219512195122 


0.858766468939860 


0.858766468941692 



By a similar argument as Lemma [III. 41 the asymptotic growth constant for the number 
of dimer- monomers on SG2,A{n) is bounded: 

31n/2,4(m)+3 lii[l+272,4M]+ln[l+372,4(m)] ^ ^ 3 In /2,4(m)+3 ln[l+2a2,4 (m)]+ln[l+3a2.4 (m.)] / , 

4x10™ — 'S'G'2,4 — 4x10™ ' V^-"/ 
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with m a positive integer. We have the following proposition. 

Proposition IV. 2 The asymptotic growth constant for the number of dimer-monomers on 
the generalized two-dimensional Sierpinski gasket SG2,i{n) in the large n limit is Zsg2a — 
0.684872262332.... 

The convergence of the upper and lower bounds is again quick. More than a hundred 
significant figures for zsg2 a can be obtained when m in Eq. fl4.6l) is equal to four. 

V. THE NUMBER OF DIMER-MONOMERS ON SGd{n) WITH d = 3,4 

In this section, we derive the asymptotic growth constants of dimer-monomers on SGd{n) 
with (i = 3,4. For the three-dimensional Sierpinski gasket SG^{n)^ we use the following 
definitions. 

Definition V.l Consider the three-dimensional Sierpinski gasket SG^{n) at stage n. (a) 
Define M^{n) = A^'^Af (5* 6*3(72)) as the number of dimer-monomers. (b) Define fsin) as the 
number of dimer-monomers such that the four outmost vertices are occupied by monomers, 
(c) Define g^ln) as the number of dimer-monomers such that one of the outmost vertices 
is occupied by a dimer and the other three outmost vertices are occupied by monomers, (d) 
Define h^{n) as the number of dimer-monomers such that two of the outmost vertices are 
occupied by monomers and the other two outmost vertices are occupied by dimers. (e) Define 
rs{n) as the number of dimer-monomers such that one of the outmost vertices is occupied by 
a monomer and the other three outmost vertices are occupied by dimers. (f) Define 53(71,) as 
the number of dimer-monomers such that all four outmost vertices are occupied by dimers. 

The quantities M3(?t,), fs^n), gs^n), h^^n), r^^n) and s^i^n) are illustrated in Fig. [9l where 
only the outmost vertices are shown. There are = 4 equivalent gs^n), = 6 equivalent 
h3{n), and (^) = 4 equivalent r3(n). By definition. 

Main) = fsin) + Ag^in) + Gh^in) + Ar^in) + s^in) . (5.1) 
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A A A A A A 



M3{n) fs{n) g^in) hsin) rs{n) 53(72) 



FIG. 9: Illustration for the spanning subgraphs M^^n), f3{n), gsin), h^(n), r^{n) and 53(71). Only the 
four outmost vertices are shown explicitly for f3{n), gs^n), h^iji), r^^n) and 53(71), where each open circle 
is occupied by a monomer and each solid circle is occupied by a dimer. 

The initial values at stage zero are /3(0) = 1, (73(0) = 0, h^i^O) = 1, r3(0) = 0, 53(0) = 3 and 
M3(0) = 10. 

The recursion relations are lengthy and given in the appendix. Some values of M^ln), 
fsin), g^in), h^^n), r^^n), s^i^n) are listed in Table IVnl These numbers grow exponentially, 
and do not have simple integer factorizations. 



TABLE VII: The first few values of Ms{n), fs{n), gs{n), /i3(n), rsin), 53(71). 



n 





1 


2 


3 


Ms{n) 


10 


945 


132,820,373,046 


49,123,375,811,021,432,878,640,796,802,876,545,882,185,505 


h{n) 


1 


51 


7,365,569,811 


2,724,928,560,954,289,860,903,291,271,266,882,549,492,483 


gsin) 





57 


7,816,070,424 


2,889,924,536,764,017,260,444,663,495,693,780,813,791,233 




1 


62 


8,289,450,499 


3,064,910,998,294,837,201,844,707,724,238,032,710,560,958 







60 


8,786,476,992 


3,250,492,861,272,219,038,243,497,885,127,347,116,333,900 


S3{n) 


3 


54 


9,307,910,577 


3,447,311,668,153,174,611,916,613,662,896,955,348,826,742 



Define 0:3 (n) = gzijn) j f^iji) and 73(72) = S3(n)/r3(?2) as in Eq. ( 13. Sp . We find {a3(n)}^]^ 
decreases monotonically and {'^z{f^)}'^=i increases monotonically with 1 < 73(71) < a3(n) for 
n>2. The values of q;3(?t,), 73(71,) and other ratios for small n are listed in Table rvTTTl 

By a similar argument as Lemma IIII.41 the asymptotic growth constant for the number 

of dimer-monomers on SG^ln) is bounded: 

ln/3(m) + 21n[l + 273(m)] In /3(m) + 2 ln[l + 2a3M] 

2x4™ - - 2x4"^ 
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TABLE VIII: The first few values of 03 (n), 73 (n) and other ratios. The last digits given are 
rounded off. 



n 


1 


2 


3 


4 




1.11764705882353 


1.06116303620220 


1.06055056935215 


1.06055052894365 


h3{n)/g3{n) 


1.08771929824561 


1.06056497054408 


1.06055052971271 


1.06055052894365 




0.96774193548387 


1.05995891923837 


1.06055049007316 


1.06055052894365 


73 (n) 


0.9 


1.05934501228135 


1.06055045043351 


1.06055052894365 



with m a positive integer. We have the following proposition. 

Proposition V.l The asymptotic growth constant for the number of dimer-monomers 
on the three-dimensional Sierpinski gasket SG^ln) in the large n limit is Zsg^ = 
0.781151467411.... 

The convergence of the upper and lower bounds is as quick as for the ordinary two 
dimensional case. More than a hundred significant figures for ^5^3 can be obtained when m 
in Eq. (15. 2p is equal to seven. 

For the four-dimensional Sierpinski gasket SGi{n), we use the following definitions. 

Definition V.2 Consider the four- dimensional Sierpinski gasket SGi{n) at stage n. (a) 
Define M^in) = N£,M{SG4^{n)) as the number of dimer-monomers. (b) Define fi{n) as the 
number of dimer-monomers such that the five outmost vertices are occupied by monomers, 
(c) Define gi{n) as the number of dimer-monomers such that one of the outmost vertices 
is occupied by a dimer and the other four outmost vertices are occupied by monomers, (d) 
Define hi{n) as the number of dimer-monomers such that two of the outmost vertices are oc- 
cupied by dimers and the other three outmost vertices are occupied by monomers, (e) Define 
ri{n) as the number of dimer-monomers such that two of the outmost vertices are occupied 
by monomers and the other three outmost vertices are occupied by dimers. (f) Define 34(77.) 
as the number of dimer-monomers such that one of the outmost vertices is occupied by a 
monomer and the other four outmost vertices are occupied by dimers. (g) Define t4^{n) as 
the number of dimer-monomers such that all five outmost vertices are occupied by dimers. 
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The quantities M4(n), /^{n), (74 (n), hii^n), r4(n), s^i^n) and t4(n) are illustrated in Fig. 
[TOl where only the outmost vertices are shown. There are (^) = 5 equivalent gi{n), (2) = 10 
equivalent h^in)^ (g) = 10 equivalent r^in)^ = 5 equivalent 54(77.). By definition, 

M^in) = fiin) + bQiin) + IQh^in) + 10r4(n) + bs^in) + ^4(77) . 

(5.3) 

The initial values at stage zero are /4(0) = 1, (74(0) = 0, h^i^) = 1, r4(0) = 0, 84(0) = 3, 
i4(0) = and M4(0) = 26. 




M4(n) /4(n) 54 (f^) /i4(n) 




r4(n) Si{n) i4(n) 



FIG. 10: Illustration for the spanning subgraphs M^ln), fiin), gi{n), h/^{n), r/^{n), 54(71) and tji{n). Only 
the five outmost vertices are shown explicitly for fiin), gi{n), h4{n), r4(n), 54(71) and ti{n), where each 
open circle is occupied by a monomer and each solid circle is occupied by a dimer. 

We write a computer program to obtain the recursion relations for SG^ln). They are too 
lengthy to be included here, and are available from the authors on request. Some values of 
M^^n), fi{n), giin), h^^n), r^^n), s^i^n), t^i^n) are listed in Table HXl These numbers grow 
exponentially, and do not have simple integer factorizations. 

Define 04 (n) = g^iji) / f^in) and 74 (n) = t^iji) / Si{n) as in Eq. (13. 8p . We find {o.iin)}'^^^ 
decreases monotonically and {'^i{n)}'^^^ increases monotonically with 1 < 74(71) < a/i{n) for 
positive integer n. The values of 04 (n), 74(71,) and other ratios for small n are listed in Table 

m 
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TABLE IX: The first few values of M4(n), /4(n), g4,{n), hi{n), r4{n), 84(11), t4(n). 



n 





1 


2 




26 


141,339 


1,567,220,397,434,550,336,692,928 


h{n) 


1 


2,460 


27,951,923,701,499,685,610,752 







3,168 


34,593,006,758,221,606,500,864 


h4,{n) 


1 


3,990 


42,806,033,106,111,666,338,688 


r4(n) 





4,852 


52,961,649,817,161,203,920,896 


S4(n) 


3 


5,683 


65,517,552,720,775,495,239,744 


t4(n) 





6,204 


81,038,847,105,336,439,783,296 



TABLE X: The first few values of 04(12), 74 (ra) and other ratios. The last digits given are rounded 
off. 



n 


1 


2 


3 


4 


04(71,) 


1.28780487804878 


1.23758948141253 


1.23734576161280 


1.23734575732423 


h4{n)/g4{n) 


1.25946969696970 


1.23741869000555 


1.23734575860203 


1.23734575732423 


r4(n)//i4(n) 


1.21604010025063 


1.23724732179398 


1.23734575559125 


1.23734575732423 


S4{n)/r4{n) 


1.17126957955482 


1.23707537334960 


1.23734575258048 


1.23734575732423 


74 (n) 


1.09167693119831 


1.23690284114717 


1.23734574956971 


1.23734575732423 



By a similar argument as Lemma lin.4l the asymptotic growth constant for the number 
of dimer-monomers on SG4{n) is bounded: 

21n/4(m) +51n[l + 274(m)] 2 ln/4(m) + 5 ln[l + 2a4M] ,n 
^^^m < , lo-4j 

with m a positive integer. We have the following proposition. 

Proposition V.2 The asymptotic growth constant for the number of dimer-monomers 
on the four- dimensional Sierpinski gasket SG4{n) in the large n limit is = 
0.876779402949.... 
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The convergence of the upper and lower bounds is as quick as for the ordinary two 
dimensional case. More than a hundred significant figures for zsga can be obtained when m 
in Eq. (15 ■4p is equal to seven. 



VI. SUMMARY 



The bounds of the asymptotic growth constants for dimer-monomers on SG2{'n), S 03(11) 
and SG^ln) given in sections lllll and IVl lead to the following conjecture for general SGd{n)- 



Conjecture VI. 1 Define ad{n) as the ratio: the number of dimer-monomers on SGd{n) 
with all but one outmost vertices covered by monomers divided by that with all outmost 
vertices covered by monomers; define 7d(n) as the ratio: the number of dimer-monomers 
on SGd{n) with all outmost vertices covered by dimers divided by that with all but one 
outmost vertices covered by dimers. The asymptotic growth constant for the number of 
dimer-monomers on the d- dimensional Sierpinski gasket SGj, is bounded 

2 In fdjm) + (rf + 1) ln[l + 27,(m)] ^ ^ 2 In /,(m) + {d + I) ln[l + 2ad{m)] 

{d + - ^^^^ - {d + 1)"^+! ■ ^ ' ^ 



We notice that the convergence of the upper and lower bounds of the asymptotic growth 
constants for dimer-monomers on SGd{n) is about the same for each integer > 2, in 



contrast to the results observed in 



35j for spanning forests on SGd{n) where the convergence 



of the bounds of the asymptotic growth constants becomes slow when d increases. 

The values of zsg^ increases as dimension d increases. Similarly for the generalized two- 
dimensional Sierpinski gasket, the values of zsG2b increases slightly as b increases. 

Compare the present results with those in Ref. 361], we find that the number of dimer- 
monomers on the Sierpinski gasket SGd{n) is less than that of spanning trees in general. 



APPENDIX A: RECURSION RELATIONS FOR 5^2,3 (n) 

We give the recursion relations for the generalized two-dimensional Sierpinski gasket 
SG2,3{n) here. For any non- negative integer n, we have 

/2,3(n+l) 
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= fl,{n) + 15fl,{n)g2,s{n) + l2fl,{n)h2,M + 84/2^,3(71)43(71) + 3/^3(71)^2,3(71) 

+mfl,{n)g2,,{n)h2M + f%{n)gl,{n) + 2Afl,{n)g2,,{n)t2M) + SS/^sH^sH 

+390/|3(7l)43(7l)/i2,3(n) + 273/13(71)^3(71) + 9/|,3(7l)/i2,3(n)t2,3(n) 
+63/|,3(n)43(7l)i2,3(n) + 180/|3(7l)^2,3(n)/l|3(^) + ^^^flMolMh^M 

+Ulf2,z{n)gl^{n) + 36/^3 (rz) (72,3 (n)/i2,3(n)t2,3(n) + 6O/I3 (71)^3 (71)^2,3 (n) 

+20/33(7l)/i3 3(n) + 264/23(7l)43(7z)/l|3(7l) + 240/2,3 (n)43(7l)/i2,3(n) + 20^3(71) 

+3/|,3(^)^2,3(^)^2,3(^) + 30/2%(7i)43(7i)/i2,3(^)*2,3(^) + ^^f2,z{n)gt,z{n)t2,z{n) 
+33/|3(7z)^2,3(n)^3H + 90/2,3(n)43(7i)/i|3(7z) + 2l43(7z)/i2,3(n) , (Al) 
g2,z{n + 1) 

= fUn)g2M + 2/|,3(7i)/i2,3(n) + 13/^3 (7z)43 (71) + flMhM 

+35/j3W52,3(n)/i2,3(n) + 60/1,3(71)43(71) + 14/^3 (7^)^2,3(r^)^2,3(r^) + 1 8/^,3 (71) /l|3(7^) 
+188/|,3(7i)43(7i)/i2,3(n) + 120/1,3(71)43(71) + llfl^{n)h2Mt2,z{n) 
+61fUn)gl,{n)t2,3{n) + 152fl^{n)g2,s{n)h^,s{n) + m fl^{n)gl^{n)h2,M 

+99/2,3(^)4,3(71) + /2',3(r^)43H + 72/|3(7l)^2,3(n)/i2,3(n)t2,3(n) 

+102/23(7i)43(7i)i2,3(n) + mUn)hl^{n) + 372/|,3(7l)43(7l)/i2 3(n) 

+310/2,3(n)43(7l)/l2,3(n) + 254,3(71) + 4/1,3(71)^2,3(^)^13(71) + 13/|,3(7l)/l|3(7l)i2,3(n) 
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+266/2,3 (n)43(n)/i|3(n) + Q^gl^{n)h2,3{n) + ^fl3{n)gl^{n)tl^{n) 
+i0fl^{n)g2,3in)hls{n)t2,3{n) + 52f2,3{n)gl-i{n)h2,3{n)t2,3{n) + 6f/|3 (71)^2,3 (^) 
+6flsin)hl,{n) + 60f2,3{n)gl,{n)hl,{n) + 3Agl,{n)hl,{n) , (A2) 

^2,3(^ + 1) 

= f2,3in)9l3in) + 443(^)y2,3(n)/i2,3(n) + 11/|,3(^)43(^) + 2/2,3(^)y2,3(n)t2,3(n) 

+4/2,3H^2,3 W + 50/1,3 W^2,3H^2,3H + W^sH^aH + 4/^^,3 (n)/l2,3Ht2,3(n) 

+21fl^{n)gl^{n)t2,3in) + 6Sfl^{n)g2,3in)hl^{n) + 191 fl^{n)gl^{n)h2,3{n) 
+56/2,3(n)43(n) + fl3in)tl3{n) + 52fl,{n)g2,3{n)h2,3{n)t2,3{n) 
+66/|,3(n)43(n)i2,3 W + 25fl,{n)hl,{n) + 289fl,{n)gl,{n)hl,{n) 
+263f2,3{n)gl,{n)h2,3{n) + 24^3(71) + 9/1,3(^)^2,3 (n)t|3(n) + 23fl,{n)hl,{n)t2,3{n) 
+'^^7fl^{n)gls{n)h2,3{n)t2,3{n) + 12f2,3{n)gt,^{n)t2,3{n) + lbQfl^{n)g2,3{n)h\^{n) 
+390/2,3(n)43(n)/i|3(n) + l^lgl^{n)h2,3{n) + ^fl3{n)h2,3{n)tl^{n) 
+19fl^{n)gl^{n)tl^{n) + 94/^3 (71)^2,3^/^2,3(^)^2,3^ + 158/2,3 (71)^3 (71) /i2,3(r^)^2,3(r^) 

+2043(7l)t2,3(n) + 20/1,3(71)^3(71) + 201/2,3(7l)^2'3(^)/i2,3(^) + 123^2',3(^)/i2,3(^) 

+Ilfl^{n)g2,3{n)h2,3{n)tl^{n) + 9/2,3(71)43(71)^^3(71) + 9/|_3(7i)/i|3(7i)t2,3(n) 
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+68/2,3(71)52,3(^)^2,3(^)^2,3(^1) + 27gl^{n)h2,3{n)t2,3{n) + 21 f2,z{n)g2,3{n)h\^^{n) 
+4l43(n)/i|3(n) , (A3) 

^2,3(^ + 1) 

= flz{n)glz{n) + 6/2%(n)43(n)/i2,3(n) + ^fl3{n)gl^{n) + 3/|3(n)43(n)t2,3(n) 
+I2fl^{n)g2,3{n)hl^{n) + hi fl^{n)gl^{n)h2,3{n) + 24/2,3(^)^3(71) 
+12/|_3(n)52,3(n)/i2,3(n)t2,3(n) + 24/|_3(n)43(n)t2,3(n) + ^flMhl^in) 
+126fl,{n)gl,{n)hl,{n) + 150/2,3(n)43(n)/i2,3(n) + 20^3(71) + 3/1,3 (71)^2,3 (n)t|3(n) 

+ 12/|,3(7l)/t2,3(7l)t2,3(^) + 99/^3 (71)^^3 (7l)/l2,3(^)^2,3(^) + 51/2,3(^)^2,3(^)^2,3(7^) 

+Ulfl,{n)g2,3{n)hl^{n) + 324/2,3 (71)^3 (71) /i|3 (71) + 120^3 (71) /i2,3(n) 
+6fl3{n)h2,s{n)tl3{n) + IS fl:,{n)gl,{n)tl^{n) + n7fl^{n)g2,3{n)hl,{n)t2,3{n) 

+ 1 86/2,3 (n)c/2,3(7l)/i2,3(7^)^2,3(7l) + 33^2,3 (71)^2,3 (t^) + SO/s 3(7l)/i2,3(^) 
+282/2,3(7l)52',3(^)^2,3(^) + 24052^,3 (7l)/i|3(7l) + /l,3(^)^2,3(^) 

+^Qfl3{n)g2,3in)h2,3in)tls{n) + 21f2,3{n)glsin)tl^{n) + 33/2%(7i)/i|3(7i)i2,3(ri) 

+ 183/2,3(^)^3(^)^3(^)^2,3(^1) + 99^2,3(^)^2,3(^^)^2,3(^1) + 87/2,3(n)52,3(ri)/i2,3(^) 

+18052',3(^)^2,3(^) + 3/13(^)^2,3(^)^13(71) + 9fl3{n)hl3inK3{n) 

+2Af2,3(n)gl^{n)h2,3(n)tl^(n) + 643(71)^^3(71) + 42/2,3 (71) 52,3 (n)/i2,3(^)*2,3(n) 
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+63gl,{n)hl,{n)t2,3{n) + 6/2,3 (n)/i^_3(n) + 3943(71)^3(71) . (A4) 
There are always 2916 = 4x3^ terms in these equations. 

APPENDIX B: RECURSION RELATIONS FOR SGsin) 

We give the recursion relations for the three-dimensional Sierpinski gasket SGs{n) here. 
For any non-negative integer n, we have 

fs{n + l) 

= fiin) + 12fi{n)g3{n) + 12fl{n)hs{n) + A8fi{n)gl{n) + Afi{n)r3{n) 
+8Afi{n)gs{n)h{n) + 72fs{n)gl{n) + 2Afl{n)gs{n)rs{n) + 30fi{n)hl{n) 
+156Mn)gl{n)h{n) + 30gt{n) + 12 fi{n)h{n)rs{n) + 36 Mn)gl{n)rs{n) 
+S4.f3{n)g3{n)hl{n) + 60gl{n)h3{n) + 24fs{n)g3{n)h3{n)r3{n) + 8gl{n)r3{n) 
+Sf3{n)hl{n) + 2igl{n)hl{n) , (Bl) 

93(n + 1) 

= fl{n)g3{n) + 3fl{n)h3{n) + 9fi{n)gl{n) + 3fl{n)r3{n) + 33fl{n)g3{n)h3{n) 

+24/3(71)53^(71) + fl{n)s3{n) + 2Afl{n)g3{n)r3{n) + 2lfl{n)hl{n) + Q6f3{n)gl{n)h3{n) 

+ I8gt{n) + Qfl{n)g3{n)s3{n) + 2lfl{n)h3{n)r3{n) + hlf3{n)gl{n)T3{n) 

+93f3{n)g3{n)hl{n) + 69gl{n)h3{n) + 3fl{n)h3{n)s3{n) + 9f3{n)gl{n)s3{n) 

+3fl{n)rl{n) + 66/3(7i)53(n)/i3(n)r3(7i) + 24y^(7i)r3(7i) + 2lf3{n)hl{n) + mgl{n)hl{n) 
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+Qf3{n)g3{n)h3{n)s3{n) + 2gl{n)s3{n) + 6f3{n)g3{n)rl{n) + 12 fs{n)hl{n)r3{n) 

+2Agl{n)h{n)rs{n) + Ugs{n)hl{n) , (B2) 
hsin + l) 

= fi{n)gl{n) + 6fi{n)gs{n)hs{n) + 6fs{n)gl{n) + 6fi{n)gs{n)rs{n) + 8/|(n)/i^(n) 
+m3{n)gi{n)hs{n) + 8gt{n) + 2fi{n)gs{n)ss{n) + Ufi{n)hs{n)rs{n) 
+S0Mn)gl{n)r3{n) + 6AMn)gs{n)hl{n) + 50gl{n)h3{n) + Afi{n)hs{n)ss{n) 
+8Mn)gl{n)ss{n) + 5/|(n)r^(n) + 80/3(71) ^3 (n)/i3(n)r3(n) + 30^f (n)r3(n) 
+26Mn)hl{n) + 87glin)hlin) + 2/|(n)r3(n)s3(n) + 16Mn)gsin)hin)ssin) 
+6gl{n)ss{n) + 18fs{n)gs{n)rl{n) + 3Afs{n)hl{n)rs{n) + 72gl{n)h3{n)r3{n) 
+Ug3{n)hl{n) + 4:fs{n)g3{n)r3{n)s3{n) + 4/3 (71)53(^1) + 8gl{n)h3{n)s3{n) 
+8f3{n)h3{n)rl{n) + 8glin)rl{n) + 2853H/i^(n)r3(n) + 4/i^(n) , (B3) 

= /sH^f (n) + 9f3{n)gl{n)h3{n) + 3gt{n) + 9f3{n)gl{n)r3{n) + 24/3(71) g3{n)hl{n) 

+27gl{n)h3{n) + 3f3{n)gl{n)s3{n) + 42/3(n)^3(n)/i3Hr3(n) + 22^|(n)r3(n) 

+18f3{n)hl{n) + 75gl{n)hl{n) + 12/3(71)^3^/13^53(71) + 6^1(71)53(71) 

+15f3{n)g3{n)rl{n) + 39f3{n)hl{n)r3{n) + 99gl{n)h3{n)r3{n) + 69g3{n)hl{n) 
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+6/3(n)^3Hr3(n)s3(n) + 9f3{n)hl{n)s3{n) + 21 gl{n)h{n) 83(71) + 21 f 3(71) h3{n)rl{n) 



+24gl{n)rl{n) + 96 g3{n)hl{n)r3{n) + 15hl{n) + 6/3(n)/i3(n)r3(n)s3(n) 

+6gl{n)r3{n)s3{n) + 12g3{n)hl{n)s3{n) + 2/3(n)rf (n) + 24g3{n)h3{n)rl{n) 

+Uhl{n)r3{n) , (B4) 
S3(n + 1) 

= gl{n) + 12gl{n)h3{7i) + 12gl{7i)r3{n) + A8gl{7i)hl{7i) + Agl{n)s3{7i) 
+84gl{n)h3{n)r3{n) + 72g3{n)hl{n) + 24^|(n) /is (n) 53(71) + 30^f (n)r^(n) 
+156g3{n)hl{n)r3{n) + 30hl{n) + 12gl{n)r3{n)s3{n) + 36g3{n)hl{n)s3{n) 

+84g3{n)h3{n)rl{n) + 60/i^(n)r3(n) + 24^3(n)/i3(n)r3(n)s3(n) + 8/i^(n)s3(n) 

+8g3{n)rl{n) + 24:hl{n)rl{n) . (B5) 
There are always 729 = 3^ terms in these equations. 
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